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Abstract. A new iterative method for non-LTE multilevel polarized radiative 
transfer in hydrogen lines is presented. Iterative methods (such as the Jacobi 
method) tend to damp out high-frequency components of the error fast, but 
converges poorly due to slow reduction of low-frequency components. The idea is 
to use a set of differently coarsed grids to reduce both the short- and long-period 
errors. This leads to the so-called multigrid (MG) methods. For the grid of N 
spatial points, the number of iterations required to solve a non-LTE transfer 
problem is of the order of 0{N). This fact could be of great importance for 
problems with fine structure and for multi-dimensional models. The efficiency of 
the so-called standard MG iteration in comparison to Jacobi iteration is shown. 
The formalism of density matrix is applied to the demonstrative example of ID, 
semi-infinite, non-magnetic, 3-principal level hydrogen atmospheric model. The 
effect of depolarizing collisions with thermal electrons is taken into account as 
well as general treatment of overlapping profiles. 



1. Formulation of the Problem 

In this paper we briefly discuss the usage of multigrid (MG ) iteration schemes to 
solve the non-LTE problem of the 2nd kind as defined by iLandi Degl'lnnocentil 
( 19871 ). Th e era of extensive development of MG methods started in 1970's by 

the work of lBrandtl ([l977). Several steps in using MG met hods applied to radia- 

five tr ansfer were made bv lSteinerl (jl99ll ) , IVathI ( 19941 ) , and lFabiani Bendicho. Trujillo Bueno. &: Auer 



( 19971 ). These authors showed that this technique leads to a great improvement 



of the convergence rate. This paper demonstrates how to apply these methods 
to a more general solution of polarized radiative transfer with realistic multi- 
level atomic models and complicated structure of overlapping lines. The effects 
of depolarizing collisions is taken into account. 

For the description of the atomic state, we adopt the density-matrix for- 
malis m and the representation in the basis of irreducible tensorial operators 
(e.g., lFano|[l957l ). The elements of atomic density matrix have the usual forrrQ 
Pg(aj), where aj is the energy level of total angular momentum j, and {k, q) are 
the multipolar components of the level {k = 0, . . . ,2j, q = —k,...,k). In sta- 



^We suppose that all coherences between different energy levels vanish due to further assump- 
tions. 
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tionary regime, the density-matrix elements are solutions of the local statistical 
equilibrium equations, 

X] ^a'j'k'q',ajkqP\{aj) = . (1) 
ajkq 

The s tr ucture of t he H- matrix has been extensively discussed by ISahal-Brechot 



([1923), i ommieil (119781 ). and Bomniicr (1980). We assume that this matrix has 



the form II = R + C, where R is the matrix of radiative rates, and C is the 
matrix of collisional rates (impact approximation). 

The radiative transfer equation for the set of four Stokes parameters S = 
(/, (5, V)^ has the usual form 

^ = J-(K-K^)5, 

ds 

where J is the emission vector of the local sources, K is the absorption matrix, 
and K'^ is the matrix of stimulated emission. All these quantities are dependent 
on radiation frequency, position vector, and direction of propagation de- 
termined by the unit vector fi. Finally, s is the parametrization of the radiation 
path along the fl direction. 



2. Standard Multigrid Method 

Most of the existing non-LTE solver s use the ni e thods based on A-operator 
splitting similar to the one of Rvbick i Hummed ( 199ll . 19921 ). Depending on 



organization, these schemes are numerically eq uivalent to the Jacobi o r Gauss- 
Seidel smoothing procedures (for details, see iPaletou &: Legerl |2005| ) . These 
smoothing procedures do reduce high frequencies of the solution fast, but poor 
convergence is achieved for low frequencies. (With "high frequencies" we mean 
those which are comparable to the spatial frequency of grid points approximating 
the continuous scale.) The principles of MG schemes are based on the idea of 
using coarse grids to reduce the low frequencies, and fine grids to smooth their 
high-frequency components. It can be showed that such a process may lead to 
the optimal CPU time demands of 0{N), N being the number of points per 
decade of optical scale. For comparison, the Jacobi and Gauss-Seidel methods 
scales approximately as 0{N'^). We have applied the non-linear version of the 
standard multig rid scheme base d on the coarse grid correction (CGC) technique 
(for details, see lHackbushlll985l ). 

CGC is the process of correction of the fine-grid approximation of the solu- 
tion using the solutions on the coarse grids. Schematically, it can be described 
in the following way: the defect (or residuum) of the fine-grid approximation is 
computed by several calls of the sweeping procedure (Jacobi, etc.); then both 
defect and the initial guess of the solution are restricted to the coarse grid, and 
a new solution on the coarse grid is obtained using these data. This coarse-grid 
solution is interpolated to the fine grid, and the density-matrix components are 
corrected. This process can be repeated recursively for every grid in order of 
increasing grid steps. This recursive process leads to the so-called V, W, or 
more complic ated diagrams, depending on the way in which the recursion is 
implemented ( Hackbushlll985l ). 
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In our code we use the parabolic short-characteristics technique of lKunasz &: Auei 
(ll98^)^and_ clS cl smoothing algorithm we use the Jacobi iteration, similar to the 
one of iManso Sainz &: Truiillo Buenol ( 20031 ). with modifications to include the 
effects of line overlapping. 



3. Detailed Model Parameters 

For demonstrating the possibilities of this technique, we chose the model of a 
semi-infinite, plane-parallel, H i atmosphere. The atmosphere is isothermal (T = 
5300 K), with constant volume density of neutral hydrogen, A'^h = lO^^cm"^, 
thermal electrons, = lO^'^'^cm"^, and protons, Np = Ng. No magnetic field is 
taken into account. The scattering of radiation is supposed to fulfill the condi- 
tions of complete frequency redistribution, and the U and V Stokes parameters 
vanish due to the symmetries of the problem. 

We adopted a hydrogen atomic model with 9 fine-structure energy levels 
(lsV2 3^^5/2) ^i^g 

non- vanishing multipole components PQ{nlj) of 
density matrix. All the coherence elements {q ^ 0) are identically zero due to 
the model's symmetries. The coherences between different energy levels have 
been neglected due to the small natural widths of the levels in comparison to 
their separations, and due to the selection rules for dipole radiative transitions. 

CoUisional rates with thermal electrons and protons for fine structure tran- 
sitions, C^/~tj^„/;/,/, were computed in part using the semiclassical theory of 
Sahal-Brechot et al. ( 19961 ). for transitions within the same shell, and in part 



using the data of collisional cross-sections from the Atomic and Molecular Data 
Information System (AMDIS; http://www-anidis.iaea.org/), for transitions 
between different shells. 

The number of logarithmically spaced nodes in the fine grid is = 257 
and 5 grids were used in total. The number of nodes in the "G — 1" grid is 
equal to Nq-i = (Nq — l)/2 -|- 1. The initial guess for the atomic density-matrix 
elements is given by the LTE populations determined in the unpolarized case. 



4. Convergence Properties and Conclusions 

Figure [1] shows a comparison between the convergence rates of the Jacobi and 
the MG methods. The effect of coarse-grid solutions is refiected by the dramatic 
evolution of the maximum relative change of the populations. The coarser the 
grid the shorter the evaluation time, and the higher the rate of approaching to 
the truncation error of the grid. The maximum relative error dominated by the 
long-period components is strongly reduced by the recursive CGC during the 
V-cycles. 

The time saving of the MG method in this particular model is about a fac- 
tor 4 compared to the Jacobi method. It must be noticed that the efficiency of 
MG increases as fast as the efficiency of the smoothing procedure. The most im- 
portant benefit from MG's is the asymptotical 0{N) behavior, which designates 
the method for use in solutions of problems that necessitate strong refin ements. 
Moreo ver t he presented ID geometry is th e slowest case as pointed out bv lSteiner 
(119911 1 and lFabiani Bendicho etHI (|l997l 1. 
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Figure 1. Convergence of the MG method with Jacobi smoothing procedure 
compared to the Jacobi method. The maximum relative change Rc of atomic 
populations, pjj, in the MG case (solid line) is compared to the Jacobi solution 
(dashed line). The maximum relative error To (with respect to the fully 
converged solution) for MG case (dotted line) and Jacobi case (dot-dashed 
line) is showed as well. The norm used is ||- • - Hoo (see lFabiani Bendicho et all 
I1997D . The graph shows the effect of 11 V-cycles with 2 pre- and 15 post- 
smoothing Jacobi iterations. 
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